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ABSTRACT 

Within a fully cosmological hydrodynamical simulation, we form a galaxy which ro- 
tates at 140 kms -1 , and is characterised by two loose spiral arms and a bar, indicative 
of a Hubble Type SBc/d galaxy. We show that our simulated galaxy has no classical 
bulge, with a pure disc profile at z=l, well after the major merging activity has ended. 
A long-lived bar subsequently forms, resulting in the formation of a secularly-formed 
"pseudo" bulge, with the final bulge-to-total light ratio B/T=0.21. We show that the 
majority of gas which loses angular momentum and falls to the central region of the 
galaxy during the merging epoch is blown back into the hot halo, with much of it 
returning later to form stars in the disc. We propose that this mechanism of redistri- 
bution of angular momentum via a galactic fountain, when coupled with the results 
from our previous study which showed why gas outflows are biased to have low angular 
momentum, can solve the angular momentum/bulgeless disc problem of the cold dark 
matter paradigm. 
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1 INTRODUCTION 

The cold dark matter paradigm successfully explains many 
properties of large scale structure observed in the Universe 
(Blumcnt hal et al.|19 84). On galactic scales, however, some 
problems exist which cast doubt on the properties of dark 
matter, and the assumption that it is non-relativistic (cold), 
not the least of which is the "angular momentum problem" . 

The angular momentum of disc galaxies is imparted by 
torques from large scale structure on dark matter and gas 
at the time at which they begin to collapse under grav- 



ity ( Peebles 1969 Barnes fc Efstathiou 19871. Disc galax- 



ies will form, under the assumption that the angular mo- 
mentum of the gas is conserved as it cools to the centre 



of the dark matter halo in which it resides (Fall & Efs 



tathiou 1980 White 1984 1. Problems with this picture of 



disc galaxy angular momentum acquisition were highlighted 
in hydrodynamical cosmological simulations, which showed 
that when gas cools to the centre of dark matter halos 
and forms stars, subsequent merging of the resultant proto- 
galaxies results in the catastrophic loss of angular momen- 



parts ( Navarro fc Steinmetz|20 00 ) . A closely related issue is 
that bulges form during mergers, as gas loses angular mo- 
mentum and rapidly cools to form stars in the centre of the 
galaxy ( Barnes &: Hernquist|199'6| . These aspects of the an- 
gular momentum problem have significantly improved with 
better modelling and resolution of feedback processes, thus 
preventing stars from rapidly forming in the central regions 



of the early collapsing dark matter halos ( Thacker & Couch- 
2001||Brook et ak]|2004| [Governato et al.||2007| ). This 



is reflected in recent simulations showing improved matches 
to the Tully- Fisher relation ( Governato et al.||2007 Agertz 
|et al.|[20T0l |Guedes et aT1|2011| see also Section |3| in this 
paper). 



Another manifestation of the problem is that the angu- 
lar momentum distribution of simulated dark matter halos 
differs markedly from that of observed galaxies. This is in 
conflict with the assumption that the distribution of the an- 
gular momentum of gas and dark matter are the same, an 
assumption which holds in non-radiative hydro-dynamical 
simulations (van den Bosch et al. 2002 Sharma & Stein- 



tum (Navarro & White 1994 1. This resulted in galaxies dom- metz|[2005 1. In particular, dark matter halos have a signif- 



inated by spheroidal components which were deficient in an- 
gular momentum when compared to their observed counter- 



icant amount of low angular momentum material ( jBullock] 
et al. 20011, and in models of galaxy formation the corre- 
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Figure 1. Face and edge on surface brightness maps in the SDSS 
g-band. Note the loose spiral arms and the existence of a bar. The 
effects of dust reprocessing have been included. 



sponding low angular momentum baryons are expected to 
form a bulge component ( |van derT Bosch 200T|. 

By contrast, many observed disc galaxies do not have 
bulge components (e.g. |Allen et al.||2006| |Kautsch||2009 



Cameron et al.|[2009 l, and have angular momentum distri- 



butions that lack the significant amounts of low angular mo- 



mentum that is found in dark matter halos (van den Bosch 



et al. 20011. The situation is worsened considerably when 
one considers that many bulges have characteristics that in- 
dicate that they were formed secularly by processes within 
the disc itself, rather than being related to merging or cold 



dark matter structure formation (Kormendy & Kennicutt 
2004} |Kormendy et al.pOlO l 



The most recent simulations of massive disc galaxies 



continue to form classical bulges, (B/T=0.19 in|Agertz et al. 


2010 and B/T=0.25 in Guedes et al.|| 


2011 see also Scan- 


napieco et al.|2010||Stinson et al. 


2010| 


Brooks et al.|201ip. 


We note that in the Guedes et al. 


(2011 


1 paper, the authors 



claim that their bulge is secular in nature, based on the low 
Sersic index of the bulge in their two component bulge+disc 
decomposition. However, they also state that their B/T has 
doubled since z = 3, indicating that their bulge is at least 
partially classical. These latest simulations have nevertheless 
formed significantly improved realisations of L, disc galax- 




Figure 2. The M s u-g colour magnitude diagram of SDSS galax- 
ies from Blanton et al. 2006, with the trend of the blue and red 
clouds indicated. Overplotted as star symbols is our simulated 
galaxy. The light blue star is the face on value, while the dark 
blue stars show various angles of inclination ( 30°, 45°, 60° and 
90° edge on, far left point). The effects of dust reprocessing are 
included, explaining the lower luminosity and redder colour of the 
simulation when viewed edge on. 



ies. Yet the angular momentum problem, and the associated 
problem of forming disc galaxies which do not have classical 
bulges, remains a serious challenge for the cold dark matter 
paradigm. 

Assuming the verity of the cold dark matter paradigm, 
it would seem that the resolution of the angular momen- 
tum problem requires either the ejection or redistribution of 
low angular momentum material. Ejection of low angular 



momentum material has certainly been proposed (Binney 



et al.||2001||Maller fc Dekel|[2002l |Dutton fc van den Bosch 



2009). Indeed, the ejection of low angular momentum gas 
by supernova explosions was shown to be key to the success 
of the first simulated bulgeless dwarf galaxies, which share 



many characteristics with observed dwarf galaxies (Gover- 
nato et al.|2010||Oh et al.|2011[ ). The processes which result 



in the ejection of low angular momentum gas from low mass 
galaxies, and the key role this has in suppressing bulge for- 
mation, were outlined in Brook et al. (2011, Paper I). 

The large scale ejection of low angular momentum gas 



meshes well with the need to enrich the IGM (Mac Low & 
Ferrara|1999}|Oppenheimer fc Dave|2006| |Shen et al.|2010[ ), 
and also with the observed low baryonic mass fraction of low 
mass galaxies (e.g. |Klypin et al. 1999| |Mandelbaum et al 



pOM IKoposov et all [20081 |Guo et al.|[2010] |Moster et al 
|2010[ > . However, in the absence of AGN feedback, more mas- 
sive galaxies with larger potential wells will have greater dif- 
ficulty ejecting their gas. This is reflected in the strength of 
the relationship between total mass and the baryonic mass 
fraction; higher mass disc galaxies have significantly higher 
baryonic mass fractions, retaining a far larger proportion of 
the cosmic baryon fraction. This means that the large scale 
outflow of low angular momentum gas is not likely to ex- 
plain the fact that many massive disc galaxies do not have 
classical bulges ( jKormendy et al.|2010 k 

Thus, processes which redistribute low angular momen- 
tum material may be necessary in higher mass disc galax- 
ies. In this paper, we propose that galactic fountains are an 
effective mechanism for such redistribution. It is generally 
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Figure 3. The rotation curve (V c = y M (< r 1 ) of the simulation 
at z = (blue line), z = 0.4 (green) and z = 1 (red). Also shown 
are the contributions from stars (dotted lines) and dark metter 
(dashed) . The rotation curve is flat at z = 1 , while it has a slight 
inner 'peak' at z = 0. 



accepted that supernova-powered bubbles drive gas out of 
the galactic disc and through the halo. The expelled gas 
eventually falls back onto the disc, a mechanism known as 
a galactic fountain (Sh apiro fc Field||1976 Bregman 19801. 
Here, we examine galactic fountains on larger scales, both 
spatially and temporally, by examining the fate of gas ex- 
pelled from the inner regions of galaxies as they are assem- 
bling their mass at high redshift. It is well established that 



supernova-driven winds can result in gas outflows ( Mathews 
fc Baker|1971 Veilleux et al.|2005 1 and that starbursts trig- 



ger strong outflows from the central regions, particularly 



during merger events (Chevalier fc Clegg||1985 



Heckman 



et al.||1990| |Strickland fc Heckman||2007| [Tremonti et ~ 
2007). Outflows are expected, and indeed observed to be 



more common at high redshift where star formation is more 



active, and potential wells shallower ( 


Madau et al. 


1996 


Pettini et al.||1998| |2000| |Simcoe et al.| 


2002; Shapley et al.| 


2003 


Adelberger et al.||2005 


I. The bipolar structure at the 


Galactic centre (Bland-Hawthorn & Cohen 2003 Su et al. 



tral regions, even locally. The fate of such outflowing gas is 
not well constrained. As explained above, there is reason to 
believe that large amounts of the gas will be lost from low 
mass galaxies, but it is likely that some or even most of the 
outflowing gas will cool back to the star forming region of 
higher mass galaxies. 

Our model is based on analysis of our fully hydro- 
dynamical cosmological galaxy formation simulations. Our 
simulation programme is running disc galaxies which span 
four orders of magnitude in stellar mass, from dwarf to L* , 
all using the same resolution and input physics. We note 
that this input physics is similar to that which simulated 



the first bulgeless dwarf galaxies (Governato et al. 20101, 
and that we reproduce all the features of those successful 
simulations in our dwarfs. The simulation presented here 
also has similar input physics and resolution to the recent 



simulation of Guedes et al. (20111, which formed a much 



improved realisation of an L* galaxy. Here we examine a 
simulated disc galaxy in an intermediate mass range, where 
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Figure 4. The Tully-Fisher relation. Open squares: data from 
Pizagno et al. 2007 using the SDSS i-band to measure Vso- Star 
symbols: our simulated galaxy, using Vgo (red star) and V m ax 
(blue star). Triangles: the recent "Eris" simulation (Guedes et al. 
2011), using V80 (red triangle) and W m ax (blue triangle). 



the majority of observed isolated galaxies are dominated by 
their disc components. 

The paper proceeds by providing details of the simu- 
lation (Section [2j|, and then describes the properties of the 
simulated galaxy (Section^, and in particular we show that 
the simulated galaxy has no classical bulge. In Section [4j we 
identify gas which is blown from the bulge region and re- 
turns via the galactic fountain to form stars at later times, 
and show that this re-accreted gas has gained angular mo- 
mentum before it forms stars. In Section \E\ we discuss our 
modelling of feedback from supernovae and its effect on the 
ISM, central energy sources and angular momentum trans- 
port. We conclude (Section [6| that galactic fountains are 
effective in redistributing the angular momentum of gas, 
and propose that this mechanism plays a crucial role in rec- 
onciling observed massive disc galaxies which do not have 
classical bulges with the ACDM paradigm. 



2 THE SIMULATION 

The simulation described here is a cosmological zoom sim- 
ulation derived from the McMaster Unbiased Galaxy Sim- 



ulations (MUG S). See Stinson et al. (2010) for a complete 



description of how those galaxies were selected and their ini- 
tial conditions were created. Here, those initial conditions 
are scaled down, so that rather than residing in a 68 ft -1 
Mpc cube, it is inside a cube with 34 h~ 1 Mpc sides. This 
resizing allows us to compare galaxies with exactly the same 
merger histories at a variety of masses. We selected a MUGs 
galaxy (15784) which had a prominent disk component for 
this study. The virial mass of the simulated galaxy in this 
study is M mr = 19.4 x 1O 1O M 

The simulation is evolved with the smoothed particle 
hydrodynamics (SPH) code GASOLINE ( jWadsley et al. 12004] ). 
GASOLINE employs cooling due to H, He, and a variety of 
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Figure 5. SDSS i-band surface brightness maps of the simulation at redshit z = 1 (left), z = 0.4 (middle) and z = (right). The disc 
grows inside out. A strong bar has developed by z = 0.4 
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Figure 6. The amplitude of the 2 nd Fourier mode of the simu- 
lation at redshift 2 = 1 (green), z = 0.4 (red) and 2 = (blue). 
The strong bar that develops by z = 0.4 is not present at z = 1. 



metal lines, and is described in detail in Shen et al. (20101 



The metal cooling grid is constructed using CLOUDY (ver- 
sion 07.02, Ferland et al. 19981, assuming ionisation equi- 



librium. A uniform ultraviolet ionising background, |Haardt| 
|fc Mada"u| ( |1996[ ) is used in order to calculate the metal cool- 
ing rates self-consistently. 

Two mechanisms are used to prevent gas from collaps- 
ing to higher densities than SPH can physically resolve: (i) 
to ensure that gas resolves the Jeans mass and does not 
artificially fragment, pressure is added to the gas accord- 



ing to Robertson & Kravtsov (20081, (ii) a maximum den- 



sity limit is imposed by setting a minimum SPH smooth- 
ing length of 0.25 times that of the gravitational softening 
length of e = 155pc. We have ~ 9 million resolution el- 
ements (TlxlO 6 gas+5.3xl0 6 dark matter+2.3xl0 6 star) 
within the virial radius at z = 0, with mean stellar particle 
mass of ~ 6400 M Q . 



2.1 Star Formation and Feedback 

When gas reaches cool temperatures (T < 15, 000 K) in a 
dense environment (n t h > 9.3 cm -3 ), it becomes eligible to 
form stars. This value for nth is the maximum density gas 
can reach using gravity, 32Af 9as /e 3 . Such gas is converted to 
stars according to the equation 



AM, 
At 



dyn 



(i) 



Here, AM, is the mass of the star particle formed, At is the 
timestep between star formation events, 0.8 Myrs in these 
simulations, M gas is the mass of the gas particle and tdyn is 
the gas particle's dynamical time, c* is the efficiency of star 
formation, in other words, the fraction of gas that will be 
converted into stars during td yn - 

Stars feed energy back into the interstellar medium 
gas where they formed. Two types of energetic feedback 
are considered in these simulations, supernovae and stel- 
lar radiation. Supernova feedback is implemented using the 



blastwave formalism described in Stinson et al. (2006) and 



deposit 10 erg of energy into the surrounding medium. 
Since this gas is dense, the energy would be quickly ra- 
diated away due to the efficient cooling. For this reason, 
cooling is disabled for particles inside the blast region 

pc and for the length of time 



R 
t - 



10' 



1.74 ciO.32 
,6.85 



-0.16 | 



-0.20 



(1977) 



it.0.32„0.34 , 



McKee & Ostriker 



04 

p o4°' 70 y r s iven in 

Here, E$\ — 10 51 erg, no is the ambient hydrogen 
4 Pofc _1 where Po is the ambient pres- 



density, and P04 = 10 
sure and k is the Boltzmann constant. Both no and Po are 
calculated using the SPH kernel for the gas particles sur- 
rounding the star. 

Metals are ejected from type II supernovae (SNII), type 
la supernovae (SNIa), and the stellar winds driven from 
asymptotic giant branch (AGB) stars. Ejected mass and 
metals are distributed to the nearest neighbour gas parti- 
cles using the smoothing kernal. Details of the manner in 
which metals are incorporated into the code are found in 



Stinson et al. (2006). The prescribed SNII and SNIa yields 



have been updated in this study. The SNII oxygen and iron 
use fits to Woosley & Weaver ( 1995 1. For SNIa, we use the 



W7 yields from Table 1 in Nomoto et al. ( 1997 1. Metal diffu- 



sion is also included, such that unresolved turbulent mixing 
is treated as a shear-dependent diffusion term, as described 
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in Shen et al. (2010). The allows proximate gas particles to 



mix their metals. Metal cooling is calculated based on the 
diffused metals. 

The effects of stellar radiation have only recently begun 
to be considered in galaxy simulations ( Hopkins et al.|20"TT I . 
Here, we model the luminosity of stars using the |Bressan| 
et al. ( 1993 1 mass-luminosity relationship. Typically, this 



relationship leads to 10 53 ergs of energy being released from 
the high mass stars over the course of the 0.8Myr between 
star formation calculations. These photons do not couple ef- 
ficiently with the surrounding ISM ( |Freyer et al. 12006] ). We 
thus do not want to couple all of this energy to the surround- 
ing gas in the simulation. To mimic this highly inefficient en- 
ergy coupling, we inject 10% of the energy as thermal energy 
in the surrounding gas, and cooling is not turned off for this 
form of energy input. It is well established that such thermal 
energy injection is highly inefficient at the spatial and tem- 
poral resolution of the type of cosmological simulations used 
here (Katz 1992 I. This is primarily due to the characteristic 
cooling timescales in the star forming regions being lower 
than the step time of the simulations. The dynamical effect 
of this feedback is negligible, as is the overall coupling of en- 
ergy to the ISM, but high mass star radiation feedback does 
reduce star formation in the region immediately surrounding 
a recently formed star particle (full details of this radiative 
feedback will be presented in a separate paper, Stinson et 
al. in prep). 



3 THE SIMULATED GALAXY: PROPERTIES 

The simulated galaxy has a total stellar mass of 1.4 x 
1O 1O M0, an HI mass of 1.9 x 1O 9 M and a total gas mass 
(within R„, r ) of 1.5x 10 10 M Q . The value log(M{H I) / L R ) = 
—0.8 fits right in the middle of the relation for observed 
galaxies Verheijen et al (2010). The disc has a metallic- 
ity which peaks at [Fe/H] — —0.5, while the cold gas has 
O/H + 12 = 8.25. The simulated galaxy's merger history is 
referred to in the analysis which follows. The interaction of 
the last major merger (mass ratio 2:1) begins at z ~ 2.7, 
with coalescence finished by z ~ 2.2 and is followed rapidly 
by 2 other significant merger events (mass ratios of 10:1 and 
20:1). We identify the period between z = 2.7 and z — 1.7 
as the galaxy's "merger epoch". Since z ~ 1.7, several dark 
matter merger events occur, but these are too low in mass 
to contain baryons. 

In Figure [T] we plot the surface brightness profile at 
z — in the SDSS g-band, shown face-on and edge-on. The 
loose spiral arms are indicative of a late type galaxy. The 
existence of a bar is also evident, and we contend that the 
simulated galaxy would be classified as SBc/d. These sur- 
face profile maps were created with Sunrise | [Jonsson||2006[ ) 
which allows us to measure the dust reprocessed spectral en- 
ergy distribution (SED) of every resolution element of our 
simulated galaxies from the far UV to the far IR with a full 
3D treatment of radiative transfer. We assume a simple re- 
lationship between metallicity of gas, which is traced in the 
simulations, and dust by using a constant dust-to metals 



fer through astronomical dust, including such effects in the 
determination of SEDs. 

The galaxy sits in the blue cloud in the M 9 , u-g colour 
magnitude diagram of SDSS galaxies, reproduced from Blan- 
(2006) in Figure [5] with the trend lines of the blue and 



ton 



red clouds indicated. The simulated galaxy is shown as star 
symbols; the light blue star shows the face-on values, while 
the dark blue stars show various angles of inclination (0°, 
30°, 45°, 60° and 90°), with greater inclination resulting in 
redder colours and lower magnitudes. The inclusion of the 
effects of dust reprocessing explains the significantly lower 
luminosity and redder colour of the simulation when seen 
edge-on. 

The rotation curve (V c = v/ -^ ) of the simulation 
is shown in Fig. [3] at z = (blue line), z = 0.4 (green) and 
z — 1 (red). Also shown are the contributions from stars 
(dotted lines) and dark matter (dashed) at those times. The 
flat curve at z = 1 reflects the success of our simulation 
in preventing excessive star formation in the central regions 
during the merger epoch. The z = curve is not perfectly 
flat, and perhaps indicates that the secular processes which 
build the bulge may be too efficient, and also indicates that, 
despite significant improvements, our simulations do not yet 
exactly reproduce the mass distribution of real galaxies. 

We plot the observed Tully-Fisher relation using data 



from Pizagno et al. (2007) who plot the SDSS i-band magni- 



tudes against Vso, the rotation velocity at the radius which 
encompass 80% of the i-band flux (open squares). Our sim- 
ulated galaxy is overplotted using star symbols, where we 
show V80 (red star) to mimic the observations, as well as 
Vmai (blue star), to show that the degree of the inner peak 
in our rotation curve is not "catasptrophic" . Also overplot- 



ted as triangles is the recent "Eris" simulation ( Guedes et al 



ratio of 0.4 Dwek (19981. Filters mimicking those on ma- 



jor telescopes are used to create mock observations. Sunrise 
uses Monte Carlo techniques to calculate radiation trans- 



20111, with colours again indicating Vso (red triangle) and 
V max (blue triangle). Using similar inputs physics and res- 
olution to our simulation, "Eris" also has a good fit to the 
Tully-Fisher relation which does not depend on the radius 
at which the measurement is made. 

The SDSS i-band surface brightness maps of the simu- 
lation at redshit z = 1 (left), z = 0.4 (middle) and z = 
(right) are shown in Fig [5] with the surface brightness range 
indicated in the colour bar. Note the inside out growth of 
the disc and the strong bar which has developed by z = 0.4. 
The 2nd mode of the Fourier harmonics of the simulations 
is shown in Fig [6] confirming that no bar exists at z = 1, 
but a strong bar with ht/hi ~ 1.2 has formed by 2 = 0.4 
(where hi is the bar length and hi is the disc scale length). 
The bar weakens somewhat by z = 0. 

In the three panels of Fig[7[ we show the decompositions 
of the SDSS i-band surface brightness maps from Fig[5] Note 
that all these one dimensional profile fits were cross-checked 
GALFIT ( |Peng et alT||2010[ ), producing similar results. At 
z = 1 (left panel), the simulation has formed a pure disc 
galaxy, best fit by a single Sersic profile having Sersic index 
n = 1.1. Recall that this is ~2 Gyrs since the last significant 
merger event. At z = 0.4 (middle panel) and z = (right 
panel) we fit a bulge+bar+disc. Also shown in the panels 
are the fit parameters, i.e. disc scale-length (h) and central 
surface brightness (u ), and the three Sersic parameters for 
the bulge and bar, i.e. the effective brightness (u e ), effective 
radius (r e ) and Sersic index (n). In our decompositions, we 
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Figure 7. Decompositions of the SDSS i-band surface brightness maps from Fig [5] At z = 1 (left panel), there is a pure disc galaxy, 
best fit by a single Sersic fit having Sersic index n = 1.1. At z = 0.4 (middle panel) and z = (right panel) we fit a bulge+bar+disc. We 
caution that the Sersic indices of the bulge components are not reliable because of gravitational softening. We do not expect that the 
derived mass associated with each component (B/T = 0,0.12,0.21 at 2 = 1,0.4,0 respectively) will be affected by this. 



derive very low Sersic indices for the bulge component. We 
caution that this quantity is not reliable, due to gravitational 
softening which will "expand" the inner mass distribution, 
with an effect out to ~ 4 softening lengths (~ 620pc in our 
simulation), which will lower the Sersic index. The amount of 
light/mass in the derived component should not be affected, 
so the derived B/T results are reliable (B/T= 0., 0.12, 0.21 
at z = 1, 0.4, respectively). 

Further, we note that other properties of secularly 
formed "pseudo" bulges are not definitive (Kormendy & 
Kennicutt 20041, but more "indicative'' 
bulge is also fiat (e bu i ge /e A 



So, although our 
ranges between 0.55 and 0.9 



depending on viewing angle in relation to the bar, where e 
is the elliptcity) and is an oblate rotator (Vmax/cv ~ 1.4, 
(-bulge = 0.48 — 0.75, again depending on viewing angle), 
by far the safest evidence that our simulated galaxy's bulge 
formed secularly is in the surface brightness profile evolution 
and its relation to merger events, from which we can deter- 
mine that the bulge did not form from gas which lost its 
angular momentum due to interactions during the merger 
epoch, but rather it formed via secular processes during the 
life of the disc. 
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Figure 8. Black line: the star formation history of our simulation. 
Blue Line: the formation history of stars which form from the 
bulge gas (BG), defined as any cold (T < 30, 000K) gas within 2 
kpc of the centre of the galaxy during the merger epoch of the 
galaxy (2.7 > z > 1.7). 20% of all stars from from BG. Redshifts 
are shown on the upper x-axis. 



4 REDISTRIBUTING ANGULAR 
MOMENTUM 

In this section, we identify and trace all gas which has 
lost angular momentum and cooled to the bulge region of 
the simulated galaxy during the galaxy's "merging epoch". 
Specifically, we identify any cold gas (T< 30, 000K) within 2 
kpc of the centre of the galaxy at any time between z = 2.7 
and z = 1.7. This gas generally leads to classical bulge for- 
mation. In what follows we refer to this gas as "bulge gas" 
(BG). There is 3.0 x 10 9 M Q of such bulge gas, which is 
~ 20% of the mass of stars in the galaxy at z = 0. Note that 
time resolution in our output steps means that it is possi- 
ble that we miss some of the cold gas which is in the bulge 
region during the merger epoch, so 3.0 x 10 9 Mq is a lower 
bound. But there is no bias in our selection; the bulge gas 
that we trace, and from which our conclusions are drawn, 
is a representative sample of cold gas in the bulge region 
during the merger epoch. 

By z = 0, 71% (2.13 x 10 9 Af Q ) of the bulge gas has 
formed stars. Of the 29% (8.7 x 1O 8 M ) remaining as gas 



at z = 0, 14% (1.2 x 1O 8 M ) has been blown entirely out 
of the galaxy, 75% (6.5 x 1O 8 M ) is in the hot halo, while 
11% (1.0 x 1O 8 M ) has subsequently cooled back down to 
form HI and is feeding star formation in the disc. In Fig|8]we 
plot the star formation history of all stars in the simulated 
galaxy at z = (black line). The blue line shows all stars 
which form from the bulge gas. Redshifts are included on 
the upper x-axis. Note the extended star formation history 
of the stars formed from bulge gas. Although 71% of bulge 
gas forms stars by z = 0, only 9% of it forms stars rapidly 
after it was first identified as being in the bulge during the 
merger epoch. The star burst from the 9% bulge gas which 
forms stars during the merger epoch ejects the remaining 
91% of the bulge gas from the bulge region, much of which 
returns to the galaxy and forms stars at a later time. We 
note that these large scale outflows also act to prevent gas 
from falling to the central region during this merger epoch, 
enhancing their effectiveness at preventing bulge formation. 

Figure [9] shows the angular momentum distribution of 
bulge gas (BG, red line) at the time that it was identified. 
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Figure 9. The angular momentum distribution of bulge gas (BG, 
red line) and stars at z = which form from the BG (blue line). 
These lines are normalised so that the area under the curves indi- 
cates the mass. Also shown is the angular momentum distribution 
of observable "baryons" at z = 0, where baryons = stars + HI 
gas (black line), normalised to have a peak at the same height as 
the bulge gas. The shape of the distribution of stars formed from 
BG (blue line) closely resembles the z = baryon distribution, 
while the BG (red line) has low angular momentum. 



The blue line shows the 2 = angular momentum distribu- 
tion of the stars which form from the BG. These two lines are 
normalised so that the area under the curves indicates the 
mass. Also shown for perspective (black line) is the angular 
momentum distribution of observable "baryons" at z = 0, 
where baryons = stars + HI gas. The baryons distribution 
is normalised to have a peak at the same hight as the bulge 
gas. The median (mean) values for the BG, stars formed 
from BG, and z = baryons are 35.7 (43.9), 92 (205) and 
164 (305) kms -1 kpc respectively.The shape of the distribu- 
tion of stars formed from BG (blue line) closely resembles 
the z = baryon distribution, while the BG (red line) has 
universally low angular momentum. The angular momentum 
of the bulge gas has been dramatically redistributed during 
the process of being blown out of the central region and cool- 
ing back to the disc via the galactic fountain. Rather than 
form bulge stars, the re-accreted gas now primarily forms 
disc stars. 

In Fig. |10| we trace the history of the 71% of the bulge 
gas which forms stars, and refer to this as "star forming 
bulge gas" (BGs). Redshifts are shown on the upper x-axis. 
In the top panel, the log of the Maximum Temperature that 
BGs reaches subsequent to the time that it was identified 
(i.e. after they were measured as being in the inner 2 kpc 
with T < 30, 000K) is plotted against the age of the stars 
which form from BGs. In the middle panel, the maximum 
distance that BGs reaches from the centre of the galaxy 
subsequent to the time that they were identified, is plotted 
against the age of BGs stars. In the bottom panel, the mean 
angular momentum of BGs stars is plotted against their age. 
9% of the BG forms stars in the bulge region shortly after 
the time that it is first identified, before it is blown out of the 
central region, and consequently forms these stars with the 
low angular momentum it had when identifed (Fig[9|. Super- 
nova feedback heats gas which has low angular momentum 
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Figure 10. Top panel: the mean of the log of the Maximum Tem- 
perature that "star forming bulge gas" (BGs, see text) reaches 
subsequent to the time that it was identified, plotted against the 
age of the stars which form from that BGs. Middle panel: the 
mean of the maximum distance that star forming bulge gas (BGs) 
reaches from the centre of the galaxy, plotted against the age of 
BGs stars. Bottom panel: the mean angular momentum of BGs 
stars is plotted against their age. Red star symbols show the me- 
dian values at each age. Supernova feedback heats gas and ejects 
it from the central region with the hottest gas travelling the fur- 
thest from the centre. The ejected gas absorbs angular momentum 
from the hot corona and later accreting gas, returning to the disk 
with a significantly different angular momentum distribution. 



and removes it from the inner central region of the galaxy, 
with the hottest gas travelling the furthest from the centre. 
The ejected gas mixes with the corona gas, as well as later 
accreting gas, and returns to the disk with a significantly 
altered angular momentum distribution. In particular, gas 
that is re-accreted from the galactic fountain at later times 
has a significanlty less negative and low angular momentum 
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baryons, a far greater amount of high angular momentum 
material, a higher median angular momentum, and forms 
stars in the disc rather than bulge. 



5 DISCUSSION 

5.1 Transition from Ejection to Redistribution 

In Paper I and in this study, we have shown that in our sim- 
ulations, there are two ways in which galaxies can alter the 
angular momentum distribution of their baryons: ejection 
in large scale outflows and redistribution of low angular mo- 
mentum material in galactic fountains. Ejection was shown 
to be a crucial process in simulated dwarf galaxies, and this 
meshes well with the low baryonic mass fractions, while re- 
distribution was shown to be a crucial process in the more 
massive simulated galaxy considered in this study. Both pro- 
cesses occur at both mass scales in our simulations. Lower 
than universal baryon fractions indicating that ejection may 
play a central role at all galaxy mass scales. 

We caution that the relative significance of ejection to 
redistribution is probably not precisely reproduced in cur- 
rent simulations. This is due to the uncertainty involving 
the inclusion of feedback and the manner in which it cou- 
ples to the ISM, a problem which is inherent to all galaxy 
formation models. In a semi-analytic study which assumes 
the ejection of low angular momentum gas, |Dutton fc van| 
|den Bosch] ( |2009[ ) found that it was necessary to couple 25% 
of SNe energy to outflows to produce realistic disc galaxy 
populations. Current estimates of the relationship between 
stellar mass and halo mass (|Mandelbaum et al.||2006| |Ko-| 
posov et al.|2008[|Guo et al.|2010||Moster et al.|2010||Dutton| 



et al.||2010||More et al.||2011[|Avila-Reese et al.||2011[ ) indi 



cate that the current simulation has a factor ~ 2 — 3 too 
many stars, implying that the role of ejection of low angular 
momentum gas was most likely underestimated for a galaxy 
of the mass considered in this study. We note that the key 
observable for detecting ejection is in fact baryon rather than 
stellar mass fraction, but the fact that our simulation has a 
reasonable value of HI/L_r means that we can assume that 
we have a similar over-abundance of baryons for our halo 
mass. What we do expect is the relative importance of red- 
stribution compared to ejection to increase in more massive 
disk galaxies, compared to their low mass counterparts, and 
that the ejection fraction due to SNe would be the lowest for 
L* galaxies where the ratio of stellar mass (Guo et al. 2010 



Moster et al. 20101 and baryonic mass (Bell et al 



2003) 



to dark halo mass are both at their highest. More massive 
systems, even groups and clusters of galaxies, also have a 



baryon fraction that is less than universal ( McCarthy et al 
20071 |Giodini et al.|2009| |Dai et al.|2(5lO| but this is widely 



attributed to the affects of Active Galactic Nuclei which are 
not included in these simulations. 



5.2 Additional Central Energy Sources 

Our simulations do not include all possible energy sources, 
and any extra central heating, for example from black holes, 



nuclear star clusters (McLaughlin et al. 2006 Rodriguez- 



|Gonzalez et al.||2011[) or a top heavy IMF in the bulge 



We again note that the outflows due to starbursts triggered 
during merger events actually prevent gas flowing to the cen- 
tral regions, enhancing their effectiveness (it was shown in 
Paper I that it is preferentially low angular momentum gas 
above and below the plane that is prevented from being ac- 
creted). 



5.3 Angular Momentum Transfer 

In a previous study of angular momentum transfer in similar 



galaxy formation simulations, Roskar et al. (20101, showed 



that outer disc warps result from torquing of the cold gas by 
a massive (misaligned) hot gas halo. Due to this torquing in 
the outer halo, the cold gas aligns with the spin of the hot 
gas halo by the time it reaches the central region. In those 
simulations, all the gas (clumpy, shocked, cold flow, recy- 
cled from disc) ended up with the same angular momentum 
distribution, due to this torquing, coupled with hydrody- 
namical drag. These processes are what is also occurring in 
our simulations. 

On the scale of local galactic fountains, significantly 
higher resolution simulations can be employed than can be 
achieved in fully cosmological studies such as ours. Hydro- 
dynamical simulations of these local galactic fountains have 
shown that both drag and accretion into clouds can ef- 
ficiently transfer angular momentum between the galactic 
corona and gas outflows (Fraternali & Binncy 2008; Mari- 
nacci et al.|2011 i. Although there are many differences be- 
tween the temperatures, times scales and velocities of the 
interactions of the galactic fountains in those studies and in 
ours, this is encouraging. Certainly, there is a need for sig- 
nificant theoretical and observation work to constrain the 
transport of angular momentum between gas with different 
origins and phases, particularly in the corona where hot halo 
gas mixes with newly accreted gas and gas blown from the 
star forming regions of the galaxy. 

We also note that disc viscosity and gas pressure within 
the disc due to supernovae and stellar winds (Valenzucla 
et al. 2007 1 may play a role in angular momentum transport, 



( Ballero et al.|2007 l, would add to the efficacy of our model. 



but some angular momentum must be imparted to the gas 
prior to re-accretion for it to fall at large enough radii to 
form part of the disc. 



6 SUMMARY 

We have formed a simulated, late type disc galaxy which has 
no classical bulge. The gas which is driven to the bulge region 
during the merger epoch, and thus has low angular momen- 
tum, is blown out by supernova powered winds which follow 
the star bursts that are triggered by the mergers. A signif- 
icant fraction of the blown out gas is re-accreted at later 
times and forms disc stars. We show that this redistribution 
of angular momentum via large scale galactic fountains can 
lead to the formation of massive disc galaxies which do not 
have classical bulges. In Paper I, we showed that galactic 
outflows will result in the ejection of low angular momen- 
tum material and hence the formation of bulgeless dwarf 
galaxies. The processes outlined in Paper I which relate to 
ejection of gas remain relevant to biases toward low angular 
momentum outflows from more massive galaxies, whether 
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or not such gas is re-accreted, and the two papers should be 
read in conjunction. 

We have shown that our self-consistent model for the 
ejection and redistribution of low angular momentum gas, 
with the latter increasing in significance in higher mass 
galaxies, can explain the existence of bulgeless disc galaxies 
that form hierarchically within a cold dark matter universe. 
Our model has the implicit assumptions that (i) energy from 
supernova during merger triggered starbursts couples effec- 
tively to the ISM and (ii) angular momentum can be trans- 
ported between between ejected gas and newly accreted and 
hot halo gas. These processes require significant theoretical 
and observational work to be better constrained, which will 
provide crucial tests of our model. 
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